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Abstract. Observations reveal nearly power-law spectta of magnetic and density 
plasma fluctuations at the subproton scales in the solar wind, which indicates the pres- 
ence of a turbulent cascade. We discuss the three-field and two-field models for micro- 
scale plasma fluctuations, and then present the results of numerical simulations of a 
two-field model of kinetic-Alfven turbulence, which models plasma motion at sub- 
proton scales. 



1. Introduction 

Magnetic turbulence is ubiquitous in astrophysical systems and it is present in labora- 
tory devices. Turbulence may be naturally generated due to various instabilities (such 
as supemovae explosions and galactic shear in the interstellar medium, tearing modes 
and shear flows in laboratory devices. Nonlinear energy cascade transfers the energy 
to smaller and smaller scales, thus distributing turbulent energy over a broad range of 
scales. At scales much larger than plasma microscales (ion cyclotron radius, skin depth, 
etc), fundamental properties of plasma turbulence can be understood in the framework 
of magnetohydrodynamics (MHD). Analytic and numerical studies of MHD turbulence 
allowed one to explain qualitatively and, in some c ases, q uantitatively in situ obser- 
vations of plasma turbuleii ce in the solar wind (e.g., iBoldyrev et aL,20ll] : IWang et al 
20Tll : IZhdankin et alll2012h . 



When the energy reaches the scales comparable to the ion Larmor radius, the char- 
acter of the turbulence changes. In laboratory flusion plasmas such micro-turbulence 
is responsible for transport phenomena. Since the large-scale guide magnetic field is 
typically strong and it cannot be easily perturbed in fusion devices, the studies have 
been mostly devoted to electrostatic fluctuations. Recently, there appeared reliable in 
situ measurements of sub-proto n plasma turbulence in the solar wind, where magnetic 
fluctuations are essential (e.g., Alexand rova et all I2009L 1201 2l : ISahraoui et all l2009l: 



Chen et al. I I2OIOL I2OI2I : ISalemetal..,2012r Such small-scale turbulence is thought to 



be responsible for energy dissipation and plasma heating in the solar wind. 

A major possibility is that significant role in subproton turbulence is played by 
kinetic-Alfven modes. Indeed, one can argue that the cascade of strong MHD tur- 
bulence (that is, turbulence of shear Alfven modes whose linearized dispersion has 
the form ta oc k^VA) is expected to transform into the cascade of kinetic-Alfven tur- 
bulence (whose linearized dispersion relation is cu oc k^kj_} at subproton scales (e.g.. 
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Schekochihin et al. I2009h . The argument goes as follows. Typical frequency of Alfvenic 
fluctuations is a> « k\\VA- The anisotropy of the energy distribution implies k\\ «; kj_ 
at small scales. At the proton gyroscale, k± ~ l/p,, one therefore estimates oj « 
kj_VA ~ kj_VTi ~ Q^i, where vji is the thermal ion velocity, Q,- is ion gyrofrequency, 
and we assumed plasma beta of order one, that is va ~ vji- Therefore, at proton gy- 
roscales, turbulence frequency is expected to be smaller than the ion gyrofrequency so 
that anisotropic kinetic Alfven modes may be effectively generated. At present, subpro- 
ton kinetic-Alfven turbulence is understood to a significantly lesser extent compared to 
its Alfvenic MHD counterpart. 

In this contribution we discuss a fluid model for micro-scale plasma fluctuations, 
and then present the results of numerical simulations of of kinetic-Alfven turbulence, 
which help to explain recent measurements of magnetic and density fluctuations at the 
subproton scales in the solar wind. 



2. Kinetic-Alfven Equations 



manv works (e.2.. Hazeltind 1983; Scott etal. 1985; 


Camarao et al. 1996 


; Terrv et al. 


2001; 


Schekochihin et al. 20091; Smith &Terrv 2011 


; Boldvrev & Pered 


20121). The 


basic assumptions are that a uniform background magnetic field is strong compared to 



magnetic fluctuations, Bq » b, and the turbulence is strongly anisotropic, k^ <«c kj_, 
where k^ and kj_ are typical wavenumbers of turbulent fluctuations in the field-parallel 
and field-perpendicular directions. To illustrate the essential physics and to set the 
notation, we start with the simplest case of small plasma beta (the ratio of thermal 
plasma energy to the magnetic energy, p = ^nn^T jB^). Then we demonstrate how the 
derived system of equations can be extended for for the case of j6 ~ 1, which is relevant 
for the solar wind stidues. 

We will assume the fluid description for the electrons. The electrons are advected 
across the magnetic field by the "E cross B" drift, v^x = cE x Bq/Bq, while their field- 
parallel motion is related to the current J\\ = -eneVe\\, and the ion parallel motion can be 
neglected. For small (3 the fluctuations of the magnetic field strength can be neglected, 
that is, the magnetic fied is represented as B = Bqz + bx- The field-perpendicular 
component is expressed through the flux function bx = z x V^, so that 7|| w Jy = 
(c/47r)Vx xbx = (c/47r)V^(/^. The flux function is the (minus) field-parallel component 
of the vector potential, ifr = -A^. 

The field-parallel force balance in the electron momentum equation gives -^\\iPe)- 
noeE\\ = 0, where the electric field is E = - - ( 1 /c)5,A. Supplementing this equation 
with the electron continuity equation, one obtains the system for the fluctuating parts 
of magnetic and density fields {J 

^l^-V + ;^^V||/., = 0, (1) 

|„^_^V0xz-Vn,-iV||7|| =0. (2) 

The electron equation can be further simplified since the electron thermal speed exceeds 
the Alfven speed, and an isothermal fluid description is possible, Tg = const. This 



'When derivatives are taken, we must distinguish the gradients along the guide field Bo, V^, from the 
gradients along the local field B = Bo + b, Vy. 
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condition applies for a collisionless plasma, and it also requires not too small plasma 
beta, p > me/mi. When collisions cannot be neglected, the electron fluid is isothermal 
if the electron diffusion time in the field-parallel direction T^iff ~ l/ik^^Vj^Tcoii) is less 
than the inverse frequencies of corresponding plasma fluctuations. 

The smallness of the plasma beta is essential for neglecting the fluctuations of 
the magnetic field strength. In the case of /J ~ 1, the fluctuations of the magnetic- 
field strength cannot be neglected, and the magnetic field is represented as B = (Bq + 
b^)z + b_L. The fact that the z-component of magnetic fluctuations should be retained 
follows from the field-perpendicular force balances in the ion and electron momentum 
equations, which can be combined to give: -VxPe-VxP(+(l/c)JxB = 0. We therefore 
derive 

- ^BoV^b, - V^pe - V^Pi - 0, (3) 

which gives an estimate b^/Bo ~ p{ne/no). For /?«:!, the fluctuations of b^ can be 
neglected, while for /3 ~ 1 they should be retained. 

It is not difficult to modify equations Q and ([T|l, taking into account b^. The 
modification comes in two ways. First, the "E cross B" velocity should be modified by 
taking into account b^, 

\e± =cEx (Bo + b,z)/B\ (4) 

where B^ k B^^ + IB^b^. Second, in the electron continuity equation Q one has to take 
into account the diamagnetic drift velocity, 

ye* = -^'^Pe^i^G + b.z). (5) 

This step requires an explanation. When the magnetic field strength does not change, 
that is, b^ - 0, it can be checked that the diamagnetic drift does not advect the electron 
density. Physically, this happens because guide centers of particles do not move when 
the diamagnetic current is present. That is why the diamagnetic drift does not enter 
Eq. ^ even in the case of a general equation of state. However, if the magnetic field 
strength changes, the magnetic curvature effects do affect the density advection, and 
terms with derivatives of ft, do not cancel out. 

Straighforward substitution of the modified drift velocities @ and ([5]l into the 
electron continuity equation then gives the modified equation ([2]): 



d_ 

di 



no Bq 



^V0xz- V 
Bo 



ne_ _ b^ 
no Bo 



eBo \nol \Bo/ eno " " 



which is derived for an arbit rary p,, but w ill be simplified using the isothermal equation 



of state (cf. iSchekochihin et al.ll2009L Eq. (C7)). In the limit of small plasma beta 



we have b^ — > 0, and Eq. Q turns into Eq. (|2]l. We should note that the field-parallel 
gradient in these equations is the gradient along the total magnetic field, that is, 

V|| - V, + -lzxV(A-V. (7) 

Bo 

In our discussion of strong kinetic Alfven turbulence, we will assume that the fluctua- 
tions are anisotropic with respect to the magnetic field in such a way that the so-called 
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critical balance between the linear and nonl inear terms is satisfied, ^nV/ ~ zx Vi/r- V; this 
condition is analog ous to k^^Bg ~ k\b (e.g.jGoldreich & Sridhad ll995UCho & Lazarian 



|2b04; Howes et al. 201 ll : iTbn Barge & Howes 2012). When this condition is satisfied, 
equations ([T] |2ll or ([1] |6]l are essentially nonlinear and three-dimensional. 

The system ([T] ^ involving the three fields, ifj, and cf) or the system ([T] |6ll 
involving the four fields rig, b^, ifj, and (p are incomplete, as they have more independent 
fields than equations. The uniqueness is restored when the systems are supplemented by 
the equations for the ions. The situation here depends on the scales considered. Above 
the ion-cyclotron scale p,- - vji/O^i, a fluid description can be justified for the ions if the 
ions are cold, which is essentially the limit of low beta. (This case is applicable for most 
laboratory plasmas, where the corresponding equations have been originally derived.) 
In this case the ions move across the magnetic field due to "E cross B" drift and the 
polarization drift, and one can write the charge concervation law, (9p/5f-i- Vj^Jj^ -i- V||7|| = 
0, where the parallel current is given by the electrons, while the perpendicular current is 
due to the polarization drift of the ions (the "E cross B" drifts are the same for ions and 
electrons, they do not lead to charge separation and do not contribute to the current). 



The resulting equation is (e.g.. lTerry et al.ll2001 ): 



— V^0-— V^xz- VV^0 
dt Bo 



- V||7||. (8) 



Equations ([T|l, Q, and (O provide the closed three-field system for evolution of electric, 
magnetic and density fields in the case of low plasma beta. 

If plasma beta is not small, the ions require kinetic description, and a simple fluid 
model is not well justified. We however will be interested in the sub-proton, disper- 
sive kinetic-Alfven waves, that is, we consider scales smaller than the ion gyroscale 
k±pi » 1. It is also convenient to introduce the ion-acoustic scale, p., = Vy/Q,- with 
Vs = {Telnii)^^'^ the ion acoustic speed. At such scales, the ions are (spatially) not 
magnetized. Moreover, we will be interested in frequencies smaller than kvTi, which 
implies the "Boltzmannian" response for the ion density fluctuations, - -ecpnolTt. 
Note that we do not require the frequencies to be smaller than the ion gyrofrequency, as 
it is implied, e.g., in gyrokinetic treatments. The quasi-neutrality condition - n^ 
then relates the electric potential to the electron density, (p = -{Tijnoe)ne. Simi- 
larly, in the three-field system, field can be removed from (|6]l according to Eq. (|3]i: 
b, = -4n{Ti + r,)«./Bo- 



3. Kinetic-Alfven Turbulence 

Let us introduce the normalized electron density and the magnetic flux function, 

n = {\+ TilT,)'l^{v,lvA)\l + {vJva)\1 + Ti/T,)]"^'^, «A - (9) 

L J no cTe 

and normalize the time and the length according to 

(i + r;/r,)i/2 



(p,/vA)[i + (v,/vA)2(i + r,vr,)] 



-^t, X = x/p,. (10) 
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We will use only the normalized variables (unless stated otherwise) and omit the over- 
tilde sign. The system ([T]), ^ then takes the form: 



dfif/ + V||W - 0, 
d,n - V||Vi«A - 0, 



(11) 
(12) 



where V|| = V^. + z x Vi/' • V^. The presented ideal system conserves the total energy E 
and the cross-correlation H, 



(13) 
(14) 



The system (II 11121 ) possesses linear waves, oc ifri^ oc exp{-iti)t + ikx). The lineariza- 
tion is done by neglecting the second term in the right-hand side of Eq. ([7]), which gives 
the dispersion relation for the kinetic- Alfven waves: 



(15) 



The linear modes are characterized by the equipartition between the density and mag- 
netic fluctuations, nk - 

For numerical simulations we supplement the equations with large-scale random 
forces that supply the energy to the system: 



dtip + V||« - TyV^ I/. + f^, 
dtn-V\^\ilj = vV\n + U 



(16) 
(17) 



The small dissipation terms serve to remove the energy at small scales (and they are 
mostly needed to stabilize the code). In a turbulent state, the energy cascades toward 
small scales while the cross-correlation cascades toward large scales. The numeri- 
cally obtained energy spectrum is shown in Fig. d] It is steeper than the spectrum 
-7/3 predicted by phenomenological theori es based on dimensional arguments (e.g., 
iBiskamp et al.|[l999l : ICho & Lazariaiil 1200 9^. It is interesting that a spectrum steeper 
than -7/3 was also inferred fror n observations of subproton magnetic and density fluc- 
tuatio ns in the solar wind (e.g., IChen et aP 120101 : lAlexandrova et al.ll2012l : IChen et aP 



I2012h . 



Various explanations have been proposed for the steeper than -7/3 spectrum of 
subproton turbulence observed in the solar wind. They include steepening of the spec- 
trum by Landau damping , weakening of tur bulence, wave-particle interactions, etc. 
(e.g.. rRudakov et al..i20ril : lHowes et al]l201ll ). In our model wave-particle interactions 
ar e absent, however, t he ste eper spectrum persists. A possible explanation proposed 



in 



Boldyrev & PerezI (120121) invoked intermittency corrections that result from two- 
dimensional structures formed by density and magnetic fluctuations. It was proposed 
that the spectrum should be close to -8/3, the value consustent with observations and 
numerical simulations. This points to an interesting possibility that the observed scaling 
is not an artifact of non-universal or dissipative effects, rather, it is an inherent property 
of the nonlinear turbulent dynamics. The spectrum may therefore be universal, analo- 
gous to the Kolmogorov spectrum of fluid turbulence. A definitive numerical study that 
requires higher numerical resolution will be conducted elsewhere. 
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Figure 1. Energy spectrum of strong kinetic -Alfven turbulence at sub-proton 
scales, obtained in two-field numerical simulations with spatial resolution 256^. 
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